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Abstract. The scattering by the inner accretion disk of 
X-ray radiation generated near the surface of a spinning 
neutron star in a low-mass X-ray binary (LMXB) has ob- 
servable effects on the waveforms of millisecond X-ray flux 
oscillations produced e.g. during type-I bursts or in the 
millisecond pulsar SAX J1808.4-3658. We study these ef- 
fects in the framework of a simplified model in which there 
is a single emitting spot on the stellar surface, which is vis- 
ible both directly and in X-rays scattered from the disk. 
The main signature of scattering from a thin disk is that 
the pulse of scattered flux leads (if the star rotates in the 
same sense as the disk) or lags (in the contrary case) the 
primary pulse of direct emission by a quarter of a spin 
cycle. This is caused by Doppler boosting of radiation 
in the sub-relativistic Keplerian flow. The disk-scattered 
flux is revealed better in energy-resolved waveforms and 
the phase dependence of the polarized flux component. 
The phenomenon discussed permits direct testing of the 
presence of standard thin disks near the neutron stars in 
LMXBs and should be observable with future X-ray tim- 
ing experiments having a few times better sensitivity than 
RXTE and also with sensitive X-ray polarimeters. 

Key words: accretion, accretion disks - polarization - 
radiative transfer - binaries: close - stars: neutron - X- 
rays: bursts 



1. Introduction 

The inner regions of the accretion disk around a weakly 
magnetic neutron star in a low-mass X-ray binary 
(LMXB) system intercept and re-emit a significant frac- 
tion of radiation generated near the surface of the neutron 
star. About half of the X-ray luminosity liberated via ac- 
cretion originates in a boundary layer of the disk or in 
a spreading layer on the stellar surface. During type-I X- 
ray bursts, the source luminosity almost entirely originates 
on the stellar surface. Therefore, scattering of the central 
emission in the inner disk should always have an impor- 
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tant effect on the properties of the observed X-ray signal 
(Lapidus & Sunyaev 1985). 

The Rossi X-ray Timing Explorer (RXTE) detections 
of millisecond X-ray flux oscillations from many sources 
have demonstrated that the neutron stars in LMXBs are 
rapidly spinning, with frequencies ranging between 300 
and 600 Hz (sec Strohmayer 2000; Van der Klis 2000 for re- 
views) . These oscillations are likely the result of azimuthal 
asymmetry in the generation of X-rays, e.g., there is a 
single emitting spot on the star (see Fig. [l]). In future X- 
ray timing experiments, a wealth of information on the 
geometrical and physical properties of neutron stars and 
the inner accretion flows in LMXBs will be obtained from 
quality measurements of oscillation waveforms. These can 
be non-sinusoidal for several reasons. First, if one consid- 
ers the situation (as is frequently done) in which there 
is no accretion disk near the neutron star, then only half 
(neglecting aberration of light) of the surface of the spin- 
ning star can be visible an any instant. Therefore, peri- 
odic eclipses of the emitting zone may take place. Second, 
since the spin velocity at the stellar equator is of the or- 
der of 0.1c (the speed of light), the emission pulse should 
be asymmetric due to Doppler boosting (the spot repeat- 
edly moves towards and away from us — see Fig. 1; Chen 
& Shaham 1989; Weinberg et al. 2000). Both general rel- 
ativistic light bending and finite spot size act to reduce 
oscillation amplitudes (Pechenick et al. 1983; Wood et al. 
1988; Riffert & Meszaros et al. 1988; Weinberg et al. 2000). 

If we take into account the accretion disk lying near 
the neutron star, it will still be possible for the above 
phase-dependent effects to take place, but additional ef- 
fects will appear at the same time. In this paper we study 
the role of the accretion disk in forming oscillation wave- 
forms. More specifically, we show that there are clear ob- 
servational effects that allow us to test the presence of a 
standard thin disk near neutron stars in LMXBs. First, 
such a disk screens half of the star from our view. Second, 
if the emitting zone is located above the disk (with respect 
to us), part of its flux is scattered by the disk in our line 
of sight. 

In the innermost regions of the disk, matter is rotat- 
ing with speeds ~ 0.5 c. Therefore, scattered photons are 
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Fig. 1. Schematic picture of scattering of radiation emitted by a spot on the surface of a spinning neutron star from 
a surrounding Keplerian accretion disk. 



Doppler-boosted along the direction of the disk rotation: 
i.e., towards us when the bright spot is irradiating the 
sector of the disk in which matter approaches us (which 
approximately corresponds to the situation presented in 
Fig. [l]) , and away from us when the opposite sector of the 
disk is being irradiated (half a cycle later). As a result, the 
phase dependence of the scattered flux component will be 
distinctly different from that of the direct emission of the 
spot. Namely, the pulse of scattered emission will lead 
the pulse of direct flux by roughly a quarter of a cycle. 
Furthermore, X-ray flux registered during this secondary 
pulse will be substantially harder and more polarized (in 
absolute units) than at other phases (Fig. ||a). These con- 
clusions only apply to binaries in which the neutron star 
rotates in the same sense as the disk. If the star rotates in 
the opposite sense (Fig. ||b) , the scattered component will 
lag the direct signal, by the same quarter of a cycle. This 
makes it possible to uncover counterrotating LMXBs (see 
Sibgatullin & Sunyaev 2000 on possibilities of formation 
of such systems). 

General relativity predicts the existence of a gap be- 
tween the last stable Keplerian orbit and the stellar sur- 
face for some (soft) neutron star equations of state (see, 
e.g., Kluzniak & Wilson 1991). In geometries with a gap, 
part of the lower stellar hemisphere is visible, and also less 
flux from the upper hemisphere is scattered in the disk. 
However, recent calculations (Sibgatullin & Sunyaev 2000) 
show that there will be no gap if the neutron star rotates 



with frequency ~ 600 Hz in the same sense as the disk. 
For the slower rotation at 300 Hz, a thin gap can exist, 
but its width proves to be smaller than the height of the 
boundary/spreading layer for the accretion rates typical of 
LMXBs (Inogamov & Sunyaev 1999; Popham & Sunyaev 
2000) . On the other hand, a large gap can arise in the case 
of counterrotation. We discuss some consequences of the 
presence of a gap in §3.2. 

A similar situation probably takes place in the mil- 
lisecond pulsar SAX J1808. 4-3658, a unique object that 
demonstrates coherent 2.5 ms flux oscillations (Wijnands 
& Van der Klis 1998; Chakrabarty & Morgan 1998; Gil- 
fanov et al. 1998). In this system, matter supposedly flows 
towards two nearly antipodal magnetic polar caps on the 
neutron star, and the accretion disk terminates at some 
distance from the stellar surface. Then, the hot spot lo- 
cated below the disk could be observable. However, the 
measured nearly sinusoidal pulse profile (Wijnands & Van 
der Klis 1998) suggests that this spot contributes very lit- 
tle, if any, to the total signal. This fact can be used for 
constraining the geometry and inclination of the system. 
The upper hot spot should be observable both directly 
and in X-rays scattered from the disk. 

Our scenario may apply to type-I X-ray bursts. Indeed, 
RXTE detections of nearly coherent oscillations during 
bursts from eight sources (see Strohmayer 2000 for a re- 
view; Wijnands et al. 2000; Galloway et al. 2000) imply 
the presence of a localized emitting zone. Apparently, this 
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Fig. 2. (a) Spin-phase dependences of various characteristics of emission from a LMXB, taking into account the effect 
of disk-scattering, for the basic model (R = R ln = 6GM/c 2 ), i = 60°, 9 = 60°, a spin frequency of / — 300 Hz, a point- 
like emitting spot (a = 0) and corotation of the star and disk. The dashed curves correspond to direct emission from 
the spot, the dash-dotted lines to the disk-scattered radiation, and the solid lines to the sum of these two components. 
The panels from top to bottom are for: photon flux, mean photon energy increment, polarized flux (in the same units 
as total flux) and polarization angle. Note that only the scattered component contributes to the polarized flux, (b) 
Same as (a), but for the case of star-disk counterrotation. 



zone slowly evolves in time (on sub-second timescales 3> 
the spin period) , as the thermonuclear burning front prop- 
agates over the stellar surface. It is not clear whether the 
innermost regions of a standard accretion disk can survive 
during a burst, when the disk is subject to strong radia- 
tion forces. In particular, radiation drag causes the matter 
flowing from the disk to lose angular momentum and spiral 



inward (Miller & Lamb 1993). Obviously, the disk has a 
greater chance to remain intact during the rise and decay 
stages of a burst, when the emitted luminosity, L, is much 
less than the critical Eddington luminosity, ^Edd- Note 
that if the inner disk is destroyed during the maximum 
of the burst, it can reconstruct very quickly at some later 
moment when the burst emission has faded substantially, 
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because the characteristic times for all processes taking 
place within several stellar radii from the star are much 
shorter than the typical burst duration (~ 10 s). 

Spots of enhanced brightness may also arise in a 
spreading layer on the surface of a neutron star (Inogamov 
& Sunyaev 1999). This layer rotates differentially, with the 
azimuthal velocity ranging between the Keplerian velocity 
(~ 0.5c) at the equator of the star, and the stellar spin 
velocity (~ 0.1c) near the edges of the layer. Under such 
circumstances, long-living whirls, analogous to typhoons, 
Jupiter's famous red spot, or "cat eyes" , may regularly ap- 
pear in two bright latitudinal belts within the layer that 
are located in the upper and lower hemispheres equidis- 
tantly from the equator. The velocity shear is largest in 
these belts. Their location is determined by the accretion 
rate, namely, they become more distant from the equa- 
tor (thus closer to the poles) with increasing accretion 
rate. Also, instabilitites developing in the accretion disk 
can modulate the accretion rate through the neck of the 
disk into the spreading layer. In particular, magnetic-field 
loops originating in the disk due to Balbus-Hawley insta- 
bility may penetrate into the layer. All these phenomena 
can lead to the formation of bright emitting zones, rotat- 
ing faster than the star itself, with quasi-frequencies up to 
the Keplerian one of ~ 2 kHz (Inogamov & Sunyaev 2000; 
Sunyaev & Revnivtsev 2000). If a spot emits a signifi- 
cant fraction of the source luminosity during several spin 
cycles, quasi-periodic flux oscillations will appear. Aver- 
aged (with a fixed period) waveforms of such incoherent 
oscillations will tend to be nearly sinusoidal. However, the 
disk-scattered flux should still manifest itself clearly in the 
power density spectrum of the total signal (i.e., in the ra- 
tio of the powers contained in different quasi-harmonics), 
as well as in the phase delays (<; 0.25) between waveforms 
measured in different energy bands and in the phase shift 
(~ 0.25) of the polarized flux component with respect to 
total flux. 

It is important to note that the considerations above 
are only applicable to atoll sources, i.e., to LMXBs with lu- 
minosities between 0.01 and 0.5 LEdd- I n this case, the ap- 
proximation of a thin accretion disk (Shakura & Sunyaev 
1973) is appropriate. At luminosities approaching LEddi 
like those of Z-sources, the boundary layer extends ra- 
dially more than one stellar radius and is geometrically 
thick (Popham & Sunyaev 2000). The adjacent disk also 
becomes thick. Therefore, our picture is no longer valid. 
Similarly in the case of L < 0.01-LEdd, Compton cooling 
becomes inefficient for matter in the disk, and a geometri- 
cally thick advection-dominated accretion flow may form 
(e.g., Narayan et al. 1998). 



thin, flat accretion disk with inner radius R m . As our basic 
model, we adopt the geometry in which R = R m = Rq = 
QGM/c 2 — 12 A km, i.e., the neutron star is the same size 
as the innermost Keplerian orbit that corresponds to the 
Schwarzschild metric. 

The spin axis is normal to the accretion plane and 
makes an angle i with the line of sight. Two types of stellar 
spin are considered: one in the same sense as the Keple- 
rian flow (corotation) and the other in the opposite sense 
(counterrotation) . 

A single homogeneously bright spot, of radius a, on 
the star emits unpolarized radiation in accordance with 
the Lambert law in its rest frame, i.e. F(fi) ~ fj, (where 
fi is the cosine of the angle between the normal to the 
spot and the viewing direction). The spot is centered at 
an angular distance 9 from the pole located in the upper 
hemisphere (see Fig. |l|). 

We also performed several computational runs, with- 
out treatment of polarization effects, for a more realistic 
situation in which the spot emits as an electron-scattering 
atmoshpere, i.e., F(n) ~ 3//(l + 2/i)/7. In this case, ra- 
diation emerges linearly polarized from the spot and is 
more beamed. We found that the effect of the additional 
beaming on the phase dependence of flux (including the 
disk-scattered component) is fairly small and likely of the 
same order as the neglected effect of gravitational light 
bending (see below) . The intrinsic polarization of the ra- 
diation emitted by the spot is also small for typical values 
of the angles i and 0, as follows directly from the classical 
results (Sobolev 1949; Chandrasekhar 1950). Indeed, the 
/^-dependence of the degree of polarization, P(n), has a 
sharp maximum of 11.7% at fi = 1. This maximum corre- 
sponds to a spot position at the limb of the star, when no 
flux is detectable. The polarized flux from the spot, i.e., 
P(/i)F(/i), peaks in the course of stellar rotation when 
\l ~ 0.5, and the maximum polarized flux, Ffi~ x , is typ- 
ically only ~ 1% of the maximum total flux, F max . The 
phase-dependent polarization effect due to disk-scattering 
proves to be much more significant (see Fig. |^), namely 
Ffnax/Ffnax ~ 5% (for the sum of direct and scattered 
radiation) for typical values of i and 6. Based on these ar- 
guments, we restrict ourselves to considering below only 
the case of a spot emitting as a Lambert source. 

Matter in the disk is in perfect circular Keplerian ro- 
tation; the velocity field (dependence on the radial coordi- 
nate r) is calculated as appropriate for the Schwarzschild 
metric: 



v<p(r) 



i?0 



-1/2 



(1) 



2. Model 

2.1. Parameters and physical assumptions 

Let a neutron star of mass M = 1.4M Q and radius R 
spin with frequency / and be surrounded by an infinitely 



so that v^Rq) — 0.5c. 

The disk is perfectly reflective, which is a reasonable 
approximation for its inner region. In fact, even if H-like 
and He-like ions of iron are abundant in the inner disk, 
these can only effectively absorb photons with energies 
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lying in a narrow band above 9 keV. For softer X-ray ra- 
diation, the reflection albedo of the inner disk should be 
almost unity. Note also that the disk's surface layer hav- 
ing a Thomson depth r J> 1 is probably completely ion- 
ized due to irradiation by X-rays from the neutron star 
(e.g., Basko et al. 1974; Nayakshin et al. 2000). For the 
above reasons, the only radiation process we consider is 
Thomson scattering. Photons incident on the disk typi- 
cally get scattered away after a few (N < 10) scatterings. 
With the temperature of matter in the disk kT ~ 1 keV 
(Shakura & Sunyaev 1973) and for incident radiation with 
E <^ 10 keV, the photon energy changes on average by a 
small fractional amount during the multiple scatterings: 
\AE/E\ ~ \N{4kT- E)/mc 2 \ < 0.1 [of course, the spec- 
trum is also Doppler-broadened: (AE/E) 2 ~ 2NkT/mc 2 ]. 
For this reason, we consider scatterings to be coherent in 
the rest frame of the Keplerian flow. To calculate the an- 
gular distribution and polarization state of scattered ra- 
diation (Thomson scattering produces linear polarization) 
as measured in the disk rest frame, we apply the classi- 
cal formulae of the theory of radiation transfer in plane- 
parallel atmospheres for the Rayleigh angular scattering 
law (Chandrasekhar 1950). 

General relativistic light bending is not taken into ac- 
count. As a result, photon trajectories are straight lines, 
which greatly simplifies the computations. At the same 
time, this simplification cannot affect the main results of 
our study qualitatively. We do not consider the redshifting 
of radiation either, because its magnitude is the same (for 
the Schwarzschild geometry) for all photons, due to their 
common origin on the stellar surface. 

All possible shadowing effects are taken into account. 
For the direct emission, these include: (1) the invisibility 
of the bright spot when it is turned away from the observer 
(provided that i + 8 > 90°), and (2) the screening by the 
disk (9 > 90°). The neutron star will, in turn, intercept 
part of the radiation scattered by the disk in the observer's 
direction if i + > 90° (the just quoted constraints are 
appropriate for the basic model with a point-like spot). 
Note that our model overestimates the duration of any 
eclipse intervals, due to the neglect of the light bending 
effect. 

2.2. The calculation procedure 

We introduce a reference frame, static with respect to the 
observer, with origin at the center of the neutron star (see 
Fig. 0). The accretion disk lies in the XY-plane, and the 
Z-axis is the spin axis of the star. For a given inclination 
angle i, the observer is at infinity in the direction 1 = 
(sini, 0, cosi). Instantaneous positions of the center of the 
emitting spot are defined by the azimuthal angle 0, in 
addition to the previously defined parameter angle 9. 

The calculation essentially consists of following, both 
in space and time, the tracks of a large number of photons 
emitted by the bright spot. Multi-dimensional integration 



is carried out over 4> (from to 2ir) , as well as over many 
zones that compose the total area of the bright spot and 
the disk. The integration over the disk is limited radially 
by an outer boundary of R ou t — 11-Ro- Only a small frac- 
tion, < 10%, of the total X-ray flux illuminating the whole 
disk is intercepted by yet farther regions (if the warping 
of the disk is neglected), causing their role in forming the 
scattered signal to be practically unimportant. It is also 
worth mentioning that for 9 > 60° (which is expected to 
occur with 1/2 probability if spot positions are random), 
more than 75% of the radiation incident on the disk is 
received by the innermost region of i?o < T < 4i?o • 

On output, we obtain the arrival time, energy shift 
(E/Eo), and polarization state for each tested photon. En- 
ergy shifts have a Doppler origin: either due to the stellar 
spin or due to the Keplerian rotation of the disk. Our code 
calculates the direct signal, the component scattered from 
the disk and the sum of these two. Several types of final 
results can be produced: (1) oscillation waveforms, i.e., 
the number of photons as a function of spin phase £ (£ = 
arrival time/spin period), (2) the phase dependence of the 
average energy increment, {E/Eq)(£) (provided that the 
spot emits monochromatic radiation of energy Eq), (3) 
oscillation waveforms in several specified narrow energy 
bands for a given energy spectrum of the spot emission, 
and (4) the phase dependence of the polarized flux com- 
ponent (namely, amplitude and the position angle of the 
polarization plane). 

The zero of £ is defined as follows. Suppose the bright 
spot is at some moment in the position (f> — and emits 
a photon in the direction 1 (along the line of sight). Imag- 
ine that at the same instant another photon is emitted 
from the center of the star, r = (0,0,0), also in the di- 
rection 1. We postulate that the second photon will ar- 
rive at the detector at phase 0. This implies that the 
first photon will reach the detector a little earlier, at 
£ = 1 — Rfc^ 1 cos(9 — i). Note that this phase difference, 
which has a time-of-flight origin, is small for / < 600 Hz: 
1 — £ < 0.025. Larger phase delays of this type can arise 
for photons that experience a disk-scattering on their way 
to the observer. The Thomson-scattering process itself is 
assumed to introduce no extra phase delay, which is true 
for our case of a geometrically thin disk. 

One phenomenon is worth mentioning here. As the 
neutron star spins, the zone of the disk that is illumi- 
nated by the bright spot rotates with a speed of 2?rfr, 
which exceeds the speed of light far enough from the star. 
For example, at / = 600 Hz, this occurs at r > 80 km. 
Such distant regions, however, intercept very little of the 
spot emission. For this reason, the contribution of super- 
luminal effects to the phase dependence of the scattered 
component is negligible. 

Let us now turn to the main physical mechanism at 
work, namely, Doppler boosting of radiation, which op- 
erates both during emission from the spinning star and 
during scattering from the disk. The disk-related effect 
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is the larger, as follows from a comparison of the Keple- 
rian velocity, ~ 0.5c, with the spin velocity at the stellar 
equator, 0.078(//300 Hz)c. Our code makes use of a set 
of equations giving the Lorentz transformations of photon 
energy E, direction of propagation k, solid angle d£L, and 
the electric (E) and magnetic (H) field vectors (Landau 
& Lifshits 1975): 



E> = E 7 ( 1 - — 



k' = 



1 



<m' = 



7(1 — kv/c) 
1 



c(7-l) ,. 
2 kv ~7 

v 



7 2 (1 — kv/c) : 



(2) 
(3) 
(4) 



Ef| = E,|, =7 [E^ 
Hj, = H y , H' x = 7 (h 

where 7 = 



[yH] 

c 

fvE" 



(5) 



-1/2 



Here, the unprimed terms are measured in the observer's 
frame, whereas the primed ones in a moving (with velocity 
v) frame, which can be either the rest frame of the bright 
spot or the rest frame of the Keplerian flow; the subscripts 
|| and _L denote a vector's components along and normal to 
v, respectively. Equations (||) govern the Lorentz transfor- 
mation of the polarization plane. Note that the polariza- 
tion degree is a Lorentz invariant (e.g., Landau & Lifshitz 
1975). 

An individual disc-scattering act is computed in three 
steps. At first, the energy, E , the direction of propaga- 
tion, ko, and the photon flux of the incident unpolarized 
radiation are transformed to the rest frame of the Ke- 
plerian flow, using Eqs. (||)-(||). Then, the photon flux 
emergent in the direction 1' and its polarization state are 
calculated in this frame, following Chandrasekhar (1950). 
Finally, the reverse transition to our frame is implemented, 
using Eqs. (H)-(||). The emergent photon energy can also 
be directly related to the incident energy: 



E = E 



1 - kpv/c 
1 - kv/c ' 



(G) 



3. Results 

3.1. Point-like emitting spot 

We begin by considering the situation in which the bright 
spot is infinitely small. Most of our results are obtained 
for the basic model (R — Ri n = Rq) and / = 300 Hz. 

In Fig. ||, we have plotted the phase dependences, 
which correspond to i = 60° and 9 = 60°, of photon flux, 



mean energy increment (for an initially monochromatic 
spectral line), polarized photon flux and polarization po- 
sition angle (which is measured relative to the meridional 
plane — defined by the normal to the disk and the line of 
sight, and takes values between —90° and 90°). The cases 
of star-disk corotation and counterrotation are compared. 

As expected, the pulse of direct emission is asymmet- 
ric, in particular, flux peaks at £ « 0.95, rather than at 
£ = 1, which would be the case if the star were spinning 
very slowly. The maximum of (E/E ) for the direct signal 
is at £ ss 0.75, which corresponds to the spot approach- 
ing us. These spin-induced Doppler effects are well-known 
(see, e.g., Weinberg et al. 2000). The spot is not visible 
directly during the eclipse centered at £ = 0.5. The direct 
emission is unpolarized by assumption. 

As regards the scattered component, flux, mean en- 
ergy increment and polarized flux all peak near £ = 0.75 
(i.e., before the peak of direct emission) in the case of 
corotation, but near £ = 0.25 (i.e., after the main peak) 
in the case of counterrotation. The partial obscuration of 
the disk by the star, which occurs around £ = 0.5, plays 
a minor role for the given values of i and 9, but becomes 
more important for larger values of these angles. It is ob- 
vious from the phase dependences for the total signal (the 
sum of the direct and scattered components) shown in 
Fig. [2] that the presence of disk-scattered flux could be 
most clearly established by means of X-ray polarimetry. 

The amplitude of the disk-scattering feature in an os- 
cillation waveform is proportional to the contribution of 
scattered flux to the time-integrated signal, which can be 
found by integrating the phase dependence. In Fig. |[ the 
relative contribution of scattered flux is shown as a func- 
tion of the inclination angle for three values of 9, and 
/ = 300. These dependences are compared with the de- 
pendences that would result if the disk were non-rotating 
and the neutron star were spinning very slowly (/ — > 0), 
so that no Doppler boosting of radiation would take place. 
One can see that the differences between the two sets of 
curves are small, i.e., Keplerian motion has a much smaller 
effect on the integral of flux than on its phase dependence. 
The spin frequency proves to be yet less important here. 
We can draw the following conclusions from Fig. |^. The 
relative contribution of the scattered component is less 
than 10% for 9 < 40°. At lower positions (with respect 
to the disk) of the bright spot, this contribution can be 
much larger (up to 100% at i = 0,9 = 90°). Provided 
that bright spot positions are random, 60° < 9 < 90° in 
50% of situations. Similarly, 60° < i < 90° for every sec- 
ond binary. Therefore, scattered emission should typically 
contribute ~ 20% to the total signal. 

Fig. |^ is similar to Fig. ||, but for the polarization de- 
gree of the phase-integrated total flux. Typically, the po- 
larization degree is ~ 3% and is somewhat enhanced owing 
to the Doppler boosting of radiation in disk-scattering. 
Our results for a non-rotating disk and slowly spinning 
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Fig. 3. Relative contribution of the disk-scattered compo- 
nent to the phase-integrated photon flux as a function of 
the angles i and 8, for the basic model with / = 300 Hz, 
a = and the case of corotation (solid lines). For compar- 
ison, the results for an imaginary non-rotating disk and 
very slowly spinning neutron star (/ = 3 Hz) are also 
shown (dashed lines). 



Fig. 4. Polarization degree of the phase-integrated total 
(direct plus scattered) flux as a function of the angles i 
and 9, for the basic model with / = 300 Hz, a = and 
the case of corotation (solid lines). For comparison, the 
results for an imaginary non-rotating disk and very slowly 
spinning neutron star (/ = 3 Hz) are also shown (dashed 
lines) . 



neutron star presented in Fig. || and Fig. [| are in agree- 
ment with the results of Lapidus & Sunyaev (1985). 

Let us now examine the effect of the scattered compo- 
nent on the harmonic content of flux oscillations. To this 
end, oscillation waveforms need to be Fourier-transformed. 
In Fig. ||, we compare the rms-amplitudes of variability 
contained in the first two harmonics. Note that much of 
the harmonic content of the direct emission is due to the 
Doppler boosting by the spinning star; more details on 
this can be found in (Weinberg et al. 2000). The scattered 
component generally causes the total signal to be more 
sinusoidal (i.e., the ratio of the amplitude of the second 
harmonic to that of the first harmonic decreases). This 
effect is mostly the result of the presence of scattered flux 
during eclipses of the bright spot (see Fig. |J) and becomes 
more pronounced for lower spot positions (see Fig. ||b). 

3.2. Effects of the disk/star size ratio, spin frequency and 
spot size 

How crucial is the condition R — R m — Rq, which de- 
fines our basic model? In Fig. ^, three different oscillation 
waveforms of disk-scattered radiation are shown, all corre- 
sponding to i = 60° and = 80° . In one case, the geometry 
is standard. In the second case, the neutron star extends 
beyond the last stable Keplerian orbit: R = R[ n > Rq. In 



first harmonic 
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Fig. 5. Amplitudes of the first and second Fourier har- 
monics in the direct flux from the spot (dashed lines) and 
in the total signal (solid lines) as a function of the inclina- 
tion angle i, for the basic model with / = 300 Hz, a = 
and the case of corotation. (a) 6 = 60°, (b) = 80°. 
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Fig. 6. Effect of the ratio of the disk inner radius, R- m , to 
the stellar radius, R, on the oscillation waveform of the 
disk-scattered component, for % = 60°, 9 — 80°, a — and 
the case of corotation. A low spin frequency, / = 3 Hz, is 
taken, so that the Doppler boosting by the star is unim- 
portant, and the direct signal (the solid line with the cor- 
responding label) is independent of R. Three geometries 
are considered: (1) the basic model (R = i?; n = 6GM/c 2 ; 
solid line), (2) R = R in = 8GM/c 2 (dotted line) and (3) 
R = 6GM/ c 2 , i? in = 8GM/c 2 (dashed line). 




Phase 



Fig. 7. Effect of the spin frequency, /, on the oscillation 
waveforms of the direct and scattered flux components, 
for the basic model with i = 60°, 9 = 80°, a = and the 
case of corotation. Two cases are considered: one of very 
slow rotation (/ = 3 Hz; solid lines), and the other of fast 
rotation (/ = 600 Hz; dashed lines). 



the third case — R = Rq < i?j n , which may approximately 
correspond to the millisecond pulsar SAX J1808. 4-3658, 
a transparent zone exists between R and i?; n . One can see 
that the amplitude of the scattered component is smaller 
by a factor of 2 in the case of this gap geometry. The wave- 
forms for the "normal-size" and "oversize" neutron stars 
are only slightly different. We conclude that the oscillation 
phase dependence and amplitude are almost independent 
of the stellar size as long as R = i?j n ;> Rq. On the other 
hand, the scattered component is sensitive to the ratio 
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Fig. 8. Effect of the bright spot size on the oscillation 
waveform of total (direct plus scattered) flux, for the basic 
model with i = 60°, 9 = 90° and / = 300 Hz. Three cases 
are compared: of corotation (short-dashed lines), of coun- 
terrotation (long-dashed lines) and also of a non-rotating 
disk (solid lines). Note the different than in the other fig- 
ures phase region, which permits a better view of the wave- 
forms. 



Ri n /R, particularly for near-equatorial spot positions, if a 
gap separates the star from the accretion disk. 

Another parameter that can have an effect on oscil- 
lation waveforms is the spin frequency, /. In Fig. ^, we 
compare two cases: one of very slow rotation, / = 3 Hz, 
and the other of fast rotation, / = 600 Hz. In the first 
case, there is no Doppler boosting of the direct emission. 
In the second case, when the equatorial spin velocity is 
0.16c, the resulting asymmetry is very pronounced, as the 
direct flux peaks at £ « 0.9. However, the scattered com- 
ponent is relatively insensitive to /. 

We have so far assumed that the bright spot is point- 
like. In reality, radiation is generated over some extended 
region of the stellar surface. Let us examine the effect of 
the spot radius, a, on oscillation waveforms. An example 
is presented in Fig. ||. A circular spot is centered at the 
equator, so that half of its area is above the disk, produc- 
ing the direct and scattered flux components, while the 
lower half is invisible and has no effect on the phase de- 
pendence. As expected, waveforms become smoother with 
increasing a. The asymmetry between the rising and de- 
clining branches of the pulse due to disk-scattering be- 
comes less pronounced as a increases. As a approaches 
90° , the waveform becomes nearly sinusoidal, and the role 
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Fig. 9. Same as Fig. g, but for the polarized flux compo- 
nent. 

of the scattered component effectively reduces to shifting 
the pulse in phase. 

Fig. H is similar to Fig. ||, but for the polarized flux 
component. From these two figures one can see that the 
pulse of polarized flux is always shifted relative to the 
pulse of total flux by approximately a quarter of a cycle, 
independent of how large the bright spot is. Whether this 
phase shift is negative or positive is determined by the 
sense of the stellar spin with respect to the disk rotation. 

3.3. Energy-resolved waveforms 

Photons arriving at the detector directly from the neutron 
star or upon scattering in the Keplerian disk have been 
Doppler-shifted in energy. The phase-dependent fractional 
shift, AE/E, is - 0.08(//300 Hz) cos [2tt(£ + 0.25)] in the 
former case (for a spot at the stellar equator and an edge- 
on view of the binary), and ~ 0.5 cos [27t(£ ± 0.25)] in the 
latter (where "+" in the argument of cos corresponds to 
the case of corotation). Therefore, energy-resolved oscil- 
lation waveforms can provide more detailed information 
on the scattering of the central emission in the disk than 
just the phase dependence of total flux. As an example, 
we consider a situation in which a spot of radius a = 60° , 
centered at the equator, emits black-body radiation with a 
temperature kT = 2 keV, a value typical for type-I bursts. 
We made computations for 5 energy bands covering the 
range E < 20 keV. The resulting oscillation waveforms are 
presented in Fig. The role of the scattered component 
apparently becomes more important with increasing en- 
ergy. Not only the relative contribution of scattered flux 




Phase 

Fig. 10. Oscillation waveforms in five energy bands for the 
basic model with i = 60°, 6 = 90°, / = 300 Hz, a = 60°, 
the case of corotation and black-body spot emission with 
kT = 2 keV. The solid lines are the total signal, and the 
dashed lines are the direct flux component. 



becomes larger, but also the pulse is much more asymmet- 
ric at high energies than at low. 

To make the analysis more quantitative, we have 
Fourier-transformed the phase dependences shown in 
Fig. [To] and thus found the phase lag of the fundamental 
harmonic as a function of photon energy (Fig. [ll]) . We see 
that the role of disk-scattering in "delaying" low-energy 
photons is more important than the corresponding role of 
stellar spin. This can be explained as follows. Both depen- 
dences plotted in Fig. [ll] approach an asymptotic value of 
—0.25 as E — > oo, because at E ^> kT flux only arises at 
phases £ ~ 0.75 (when both the spot and the illumination 
front on the disk are moving towards us). The phase lead 
will constitute a significant fraction of this 0.25 value for 
E £ kT/(AE/E) [because of the exp (-E/kT) factor in 
the spectral law]. Therefore, in order to have a substantial 
hard-to-soft phase difference due to disk-scattering, one 
needs E £ 2fcT, as compared to E > 10(300 Hz//)fcT 
in the case of stellar rotation. Thus, in the latter case, 
large phase leads (~ 0.1) can only materialize very far in 
the Wien tail of a black-body spectrum, i.e., in a region 
where virtually no emission is generated. Note that in the 
case of a counterrotating disk, hard emission will lag soft 
emission. 
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Fig. 11. Phase lag of the fundamental harmonic as a func- 
tion of photon energy derived from the oscillation wave- 
forms shown in Fig. [h]. The solid bars correspond to the 
total signal, and the dashed bars to the direct flux com- 
ponent. 



4. Conclusion 

We have considered the effect of X-ray scattering by a thin 
accretion disk on the phase dependence of spin-modulated 
emission from a neutron star with a weak magnetic field. 
This effect makes it possible to test the presence of stan- 
dard thin accretion disks around the neutron stars in 
LMXBs. The only other existing experimental method 
to probe Keplerian accretion disks near compact objects 
(particularly black holes) is through studying the profiles 
and variability of iron fluorescence lines (see Fabian et 
al. 2000; Reynolds 2000 for recent reviews). Furthermore, 
the new method provides a unique opportunity to check 
whether the neutron star is corotating with the disk, and 
potentially even to determine the misalignment angle be- 
tween the spin axis and the disk symmetry axis. Such 
measurements can significantly enrich our knowledge of 
the origin and evolution of LMXBs. 

The results of this paper demonstrate that the disk- 
scattered flux component manifests itself better in energy- 
resolved oscillation waveforms. Nearly coherent oscilla- 
tions during type-I bursts, especially those detected dur- 
ing the burst rise, are probably the best candidates to look 
for signatures of disk-scattering. It is possible that RXTE 
data can already be used for such a search. For example, 
the phase lag of low-energy photons during a burst from 
Aql X-l reported by Ford (1999) may be the result of 
Doppler boosting in a thin disk. It is important that this 



interpretation is self-consistent, as opposed to the alterna- 
tive explanation based on Doppler boosting by the neutron 
star, as shown in §3.3. The case of the millisecond pulsar 
SAX J1808. 4-3658, for which similar phase lags and also 
distortions of the pulse profile have been detected (Cui et 
al. 1998; Revnivtsev 1999; Ford 2000), appears to be more 
complicated. The principal difficulty here is that the ob- 
served energy spectrum has an approximately power-law 
shape. If this spectral distribution originates at the hot 
spot, then no energy dependence of the pulse profile can 
arise as a result of Doppler boosting (see, e.g., Chen & 
Shaham 1989). 

More detailed studies of the waveforms of X-ray flux 
oscillation in LMXBs could be possible with ~ 6 m 2 X-ray 
detectors (i.e., 10 times larger than RXTE), such as the re- 
cently proposed EXTRA and RAE missions (Barret 2000; 
Kaaret 2000; http://www.cesr.fr/~barret/extra.htm]). 



A particularly promising way to look for spin-phase 
signatures of disk-scattering, at least in the millisecond 
pulsar SAX J1808. 4-3658, is by means of sensitive 
X-ray polarimetry. The case of SAX J1808.4-3658 (or 
similar objects which may be discovered in the future) 
is favourable because this source demonstrates regular 
pulsations during long time intervals, which allows one to 
accumulate a statistically large amount of data. On the 
contrary, such statistics is hardly achievable for nearly 
coherent burst oscillations, or for possible quasi-periodic 
oscillations originating in the spreading layer. Previously, 
only the time-integrated polarization due to Thomson 
scattering in the accretion disk (without allowance for 
its rotation) was discussed as an observational task for 
near-future X-ray polarimeters (Meszaros et al. 1988). 
Such polarimeter will be part of the scientific payload of 
the Spectrum- X-Gamma observatory (Kaaret et al. 1993; 



http : / /hea- www . harvard . edu / ~ kaaret /sxrp /index . html ) . 
Also, NASA is currently attaching great importance to 
a programme of long duration balloon flights, lasting up 
to 100 days (Jones 2000). This programme will no doubt 
rekindle interest to balloon flights of X-ray polarimeters. 
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